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Abstract We review recent computational results for hexagon patterns in non- 
Boussinesq convection. For sufficiently strong dependence of the fluid parameters 
on the temperature we find reentrance of steady hexagons, i.e. while near onset the 
hexagon patterns become unstable to rolls as usually, they become again stable 
in the strongly nonlinear regime. If the convection apparatus is rotated about a 
vertical axis the transition from hexagons to rolls is replaced by a Hopf bifurcation 
to whirling hexagons. For weak non-Boussinesq effects they display defect chaos 
of the type described by the two-dimensional complex Ginzburg-Landau equa- 
tion. For stronger non-Boussinesq effects the Hopf bifurcation becomes subcritical 
and localized bursting of the whirling amplitude is found. In this regime the cou- 
| pling of the whirling amplitude to (small) deformations of the hexagon lattice 

becomes important. For yet stronger non-Boussinesq effects this coupling breaks 
C^*) • up the hexagon lattice and strongly disordered states characterized by whirling 

£f) ' and lattice defects are obtained. 

o ' 

O " 1 Introduction 



Rayleigh-Benard convection has been a paradigmatic system for studies of spontaneously form- 
ing spatial or spatio-temporal patterns (e.g. pQ) far from equilibrium. In recent years exciting 
results have been obtained for the stability and dynamics of structures that are connected with 
roll convection. In particular, the spiral-defect chaos obtained in convection of fluids with low 
Prandtl number [2 3 4|5|6|7] and domain chaos |8|9|1Q|11] driven by the Kiippers-Lortz insta- 



bility [12] in rotating systems have attracted notable interest. In most of these investigations 
the systems have been kept close to the Oberbeck-Boussinesq approximation by minimizing the 
dependence of the fluid parameters on the temperature to avoid the appearance of cellular or 
hexagonal structures. 

Variations of the fluid parameters with the temperature, i.e. non-Oberbeck-Boussinesq ef- 
fects, break the up-down symmetry and give rise to hexagons at the onset of the convection. 
Weakly nonlinear analysis shows that typically hexagons become unstable to rolls further above 
threshold [13J. This scenario has been confirmed in quite a number of experimental investiga- 
tions [lj and it has often been assumed that in convection hexagonal patterns are confined to 
the regime close to onset. Presumably due to this assumption, hexagonal convection in strongly 
nonlinear non-Boussinesq convection has not received as much attention. 

Recently a couple of experiments |14)15j showed that even in the strongly nonlinear regime 
hexagonal convection patterns can be observed stably. In one case hexagons were observed at rel- 
atively high Rayleigh numbers (e = (R—R c )/R c « 3.5) under conditions in which the Oberbeck- 
Boussinesq approximation was quite well satisfied. Correspondingly, hexagons with up- and with 
down- flow in the center were observed to coexist in adjacent domains [14]. Subsequently, the 
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linear stability of up- and down- hexagons in this regime of Boussinesq-convection was confirmed 
by a numerical stability analysis |16|17j . In the other case a strongly non-Boussinesq situation 
was investigated and hexagons were observed not only close to threshold but also again as 
reentrant hexagons at higher Rayleigh numbers, e = 0(1) [15J. As the non-Boussinesq effects 
were increased the intermediate e-range over which rolls were the preferred planform shrank 
and eventually hexagons were found to dominate rolls from directly at onset all the way to 
e = (9(1). The restabilization of hexagons at larger e was attributed to the high compressibility 
of the fluid; the experiment was performed with SF6 near its thermodynamical critical point 

If the chiral symmetry of the convective apparatus is broken by rotating it about the ver- 
tical axis the Kiippers-Lortz instability [12] leads to new complex dynamics. Above a critical 
frequency rate domain chaos sets in, which is characterized by a persistent switching between 
roll patches with different orientation. In the non-Boussinesq case hexagons are steady near 
threshold. However weakly non-linear theory predicts a secondary Hopf bifurcation to oscilla- 
tions in which the three modes forming the hexagonal lattice oscillate with the same frequency 
but with a phase shift of 2tt/3 |18J19| . 

Here we review recent computational results for strongly nonlinear non-Boussinesq con- 
vection 20 21 22 and illustrate them further with additional data. The paper is organized as 
follows. After outlining the basic equations for non-Boussinesq convection (Section [2]) we dis- 
cuss reentrant hexagons in fluids with high and to some extent also with low Prandtl number 
(section [3}. The reentrance is clearly illustrated with numerical simulations in relatively large 
systems. In Section |4] we address non-Boussinesq convection in rotating systems. We review 
the transition to whirling hexagons in the super- and the subcritical case, present new results 
for a period-doubled state, and illustrate the mechanism underlying the apparently persistent 
whirling chaos, which involves defects in the hexagon lattice. Conclusions are given in Section 

El 



2 Basic equations 



The basic equations that we use for the description of non-Boussinesq convection have been 
discussed in detail previously |23|21j . We therefore sketch here only a brief summary. We con- 
sider a horizontal fluid layer of thickness ft, density p, kinematic viscosity heat conductivity 
A, thermal diffusivity and specific heat c p . The system is heated from below (at temperature 
T\) and cooled from above (at temperature T2 < X\), and in the rotating case it rotates about 
a vertical axis with angular velocity otherwise uo = 0. 

To render the governing equations and boundary conditions dimensionless we choose the 
length ft, the time ft 2 /^o 7 the velocity Ko/d, the pressure po^o^o/^ 2 , and the temperature 
T s = voKo/aogh 3 as the respective scales. The subscripted quantities refer to the respective 
values at half depth of the fluid layer in the conductive state. The non-dimensionalization gives 
rise to three dimensionless quantities: the Prandtl number Pr = ^o/^o, the Rayleigh number 
R = ao ATgh 3 /voKo, and the dimensionless rotation rate Q = ouh 2 /vQ. Furthermore, we write 
the equations in terms of the dimensionless momentum density V{ = phiii/po^o instead of the 
velocities U{. The dimensionless form of the temperature T = T/T S: heat conductivity A = A/Ao, 
density p = p/po, kinematic viscosity v = v/vq, and specific heat c p = c p /c p o will be used in 
the ensuing equations and the hats dropped for simplicity. In dimensionless form the equations 
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for the momentum, mass conservation and heat are then given, respectively, by 



= -dip 



P 



with the dimensionless boundary conditions 
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e. 



+S i3 [ l + li(-2z + -) ) 



"H*(t)H(t) 



+2f2e ij3 Vj 



diVj = 0, 



1 

pc p 



R^(l+ l3 z). 
P 



(2) 



(3) 



v(x, y, z, t) = 0(x, y, z, t) = at z = ±-. 



(4) 



Here O is the deviation of the temperature field from the basic conductive profile, Sij is the 
Kronecker delta, and e^3 the Levi-Civita tensor. Summation over repeated indices is assumed. 

In the rotating case, we assume small rotation rates and neglect the centrifugal force term 
relative to gravity, i.e., we assume i? 2 r « g. Notice that rotation breaks the chiral symmetry 
through the presence of the Coriolis term 217 x u. 

We consider the non-Boussinesq effects to be weak and retain in a Taylor expansion of all 
material properties only the leading-order temperature dependence beyond the OB approxi- 
mation. For the density this implies also a quadratic term with coefficient 71. It contributes, 
however, only to the buoyancy term in (pQ); in all other expressions it would constitute only a 
quadratic correction to the leading-order non-Boussinesq effect. Thus, the remaining tempera- 
ture dependence of the fluid parameters p, ^, A, and c p in (|l|2|3p is taken to be linear 



p(©) = l-7o(-*+|), 
K0) = l+72(-^+|), 
A(0) = l + 73 (-^+|), 
c p (0) =l + 74 (-z+|). 



(5) 
(6) 

(7) 

(8) 



The coefficients 7^ give the difference of the respective fluid properties across the layer. They 
depend therefore linearly on the Rayleigh number, 



7i (AT) = 7? 



7 c (l + e), 



(9) 



where 7? is the value of 7^ at the onset of convection and e = {R— R c {p(^)) / R c {p(f) is the reduced 
Rayleigh number. 

The linear analysis yields the critical Rayleigh number R c as well as the critical wavenumber 
q c . Both depend on the non-Boussinesq-coefficients jf which in turn depend on R c . Thus, in 
principle one obtains an implicit equation for the 7 2 c . The shift in the critical Rayleigh number 
away from the classical value R c = 1708 due to the non-Boussinesq-effects is, however, quite 
small (less than 1 percent) and therefore the resulting change in the 7? is negligible. In this 
paper we therefore choose the 7 2 c corresponding to R c — 1708. 
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Figure 1. (Color online). Phase diagram showing the results of simulations for a layer of water with 
thickness h — 1.8 mm in a circular cell with diameter L = 16 x 2n/q c . The simulations have been 
carried out for To = 24 °C and 28 °C (for further parameters see Figs i2l3|) . Circles correspond to rolls 
and triangles to hexagons. The full and dashed line correspond to the stability limits for amplitude 
instabilities of hexagons and rolls, respectively, as determined by the Galerkin calculation. 



Details of the numerical scheme used to determine the stability of patterns in the strongly 
non-linear regime, and for the integration of the Navier-Stokes equations can be found in ref- 
erence [21]. 



3 Reentrant hexagons 

The classic weakly nonlinear analysis for weakly non-Boussinesq convection [13J predicts that 
not far above threshold the hexagons become unstable to a mixed-mode, which in turn is also 
unstable, resulting in a transition to roll convection. This transition has been addressed in ex- 
perimental studies in CO2 [24] and water [26 [25], respectively. In both studies a transition to 
rolls was observed. The quantitative agreement with the weakly nonlinear theory was, how- 
ever, limited. Among the reasons identified for the poor agreement were the boundaries, which 
enhance the tendency towards rolls. 

In recent work on the stability of fully nonlinear hexagon solution of the Navier-Stokes 
equations we have again addressed the transition from hexagons to rolls [21 20J and found that 
a further cause for disagreement between weakly nonlinear theory and experiments may be 
given by the fact that even for seemingly moderate non-Boussinesq effects the transition from 
hexagons to rolls can be shifted towards larger values of e and can, in fact, disappear altogether. 

Results of our stability calculation for a layer of water of height h = 1.8 mm are shown as a 
function of the mean temperature To of the layer in FigHJ The mean temperature is denned as 
To = (Ti + T2)/2, where T\ and T2 are the temperatures at the bottom and top plates respec- 
tively, and To is kept constant in the simulations of the Navier-Stokes and energy equations. 
Hexagons are stable below and to the left of the solid line, whereas rolls are linearly stable to 
hexagons above the dashed line. With decreasing mean temperature (and with decreasing layer 
height) the critical temperature for the onset of convection increases and with it the variation 
of the fluid parameters across the layer. In addition, for lower mean temperatures the temper- 
ature dependence of the density of water develops a strong quadratic contribution. Thus, with 
decreasing mean temperature the non-Boussinesq effects increase. This shifts the line indicating 
the instability of hexagons to rolls to larger values of e. More importantly yet, we found that 
the hexagons can regain stability with respect to rolls if the Rayleigh number is increased well 
beyond the weakly nonlinear regime. This line of reentrance is shifted towards lower value of 
e as the non-Boussinesq effects become more pronounced and in fact merges eventually with 
the lower stability limit of the hexagons at a temperature T m . Below T m the hexagons do not 
undergo any instability with respect to rolls over the whole range of e shown. 

More detailed stability analyses of the hexagons with periodic boundary conditions using a 
general Galerkin-Floquet ansatz to capture side-band instabilities show that the only relevant 
side-band perturbations are long- wave and steady, as is the case in the weakly nonlinear regime. 
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Figure 2. Series of snapshots in a circular cell of water of thickness h — 1.8 mm and mean temperature 
To = 24 °C. The diameter of the cell is L = 16 x 27v/q c , and the snapshots correspond to an integration 
time of 600 t v . The corresponding non-Boussinesq coefficients are Jq = 0.0039, ji = 0.3326, 7§ = 
-0.3641, 7§ = 0.0467, 7I = -0.0026, and the Prandtl number Pr = 6.3. 



The range of stable wavenumbers q shifts towards smaller values of q as e is increased. This 
reflects the fact that even in the absence of non-Boussinesq effects hexagons can become linearly 
stable relatively far from threshold, e > 2. These hexagons are, however, stable with respect to 
side-band instabilities only for relatively low values of q [T6j . 

A seemingly similar reentrance of hexagon convection has been observed by Roy and Stein- 
berg in SF6 near the thermodynamical critical point [15]. There it was argued that since the 
non-Boussinesq effects in that system are not very large near onset the reentrance is due to 
the large compressibility of the fluid in this parameter regime. By assuming that the working 
fluid is incompressible, which is an excellent approximation for water, our computations show 
that high compressibility is not needed for reentrance. The significance of the compressibility 
for the occurrence of the reentrance has also been called into question in [27 28J. Moreover, 
the wavenumber of the reentrant hexagons of [15J does not decrease with increasing e, but 
stays quite close to q C: whereas a characteristic feature of the reentrant hexagons found in our 
computation is that the wavenumber is significantly lower than q c (cf. [16]). 

With respect to the experiments in water [25] our stability computations show that even in 
that system reentrance should have been accessible just slightly above the Rayleigh numbers 
investigated in those experiments (0 < e < 0.14). For the experiment in CO2 [24] we find that 
the parameters correspond to the regime in which hexagons do not become unstable to rolls at 
all up to e = 1, similar to the situation shown in Fig Q] for water for To < T m . 

To address the influence of the side walls we have performed numerical simulations of the 
Navier- Stokes equations in which we mimic the circular sidewalls by applying a strong radial 
subcritical ramp in the Rayleigh number that suppress any convection outside a certain radius 
[3]. The system size is chosen to be twice as large as in [20 , i.e. here L = 16 x 2tt jq c . Random 
initial conditions are used. Fig. [2] shows snapshots of simulations for a mean temperature of 
To = 24° C. Consistent with the phase diagram shown in Fig Q] no transition from hexagons to 
rolls is found. As expected from the shift of the stable wavenumber band towards smaller value 
of q the hexagon wavenumber is decreasing with increasing e. Note that the defects that cause 
the change in wavenumber do not trigger a transition to rolls. 



For larger mean temperature, To = 28°C, the non-Boussinesq effects are weaker. As a 
consequence, for intermediate values of e rolls invade the hexagon pattern from the side walls 
similar to the experiments of [25 24J and for values of e in the range 0.6 < e < 0.1 disordered 
roll patterns are obtained, which are characteristically perpendicular to the sidewalls (Fig|3]). 
For e = 1.2 the hexagons take over again and we find a new interesting scenario: coexistence 
between rolls and reentrant hexagons. This state was not obtained in our previous simulations 
of water due to the small aspect ratio used [21 J nor in simulations of reaction diffusion system 
[29] . At e = 1.4 the non-Boussinesq effects become strong enough to remove all the rolls and 
the reentrant hexagons take up the whole convective cell. 

The results of the simulations are summarized in FigQ] with triangles denoting hexagons 
and circles marking roll states. Compared to Fig. 9 in [21] the regime of disordered roll patterns 
is shifted towards larger To, i.e towards weaker non- Boussinesq effects. This is due to the larger 



6 



Will be inserted by the editor 




£=0.02 £=0.1 £=0.6 




£=1.0 £=1.2 e=1.4 

Figure 3. Sequence of snapshots for To = 28 °C in a circular cell of water of thickness h — 1.8 mm. 
The diameter of the cell is L = 16 x 27i/q c , and the snapshots correspond to an integration time of 
300 t v . The corresponding non-Boussinesq coefficients are 7§ = 0.0036, 71 = 0.2122, 7§ = —0.2725, 
7 § = 0.0352, 7! = -0.0013, and the Prandtl number Pr = 5.7. 

convection cell used in the present simulations. This reduces the influence of the side walls. They 
enhance the nucleation of rolls, which then can invade the hexagonal pattern in the bulk. 



4 Convection with rotation 

If the whole system is rotated about a vertical axis the Coriolis force breaks the chiral sym- 
metry, which transforms the steady bifurcation off the hexagons to the unstable mixed mode 
into a Hopf bifurcation to whirling hexagons [18,19]. In [3Q|31] the complex Ginzburg-Landau 
equation describing spatial modulations of the complex oscillation amplitude TL of the whirling 
hexagons has been derived from three coupled Ginzburg-Landau equations that provide the 
weakly nonlinear description of the hexagons, 

d T H = aH + dV 2 H - cH\H\ 2 . (10) 

In principle, this CGL is still coupled to the phase modes of the underlying hexagon pattern. 
However, at the band center q c the coupling coefficients vanish, irrespective of the coefficients of 
the three coupled Ginzburg-Landau equations. Moreover, the imaginary parts of the coefficients 
d and c are always in the regime of bistability of the homogeneous oscillations and defect chaos 

Recently, we have investigated this Hopf bifurcation within the fully nonlinear Navier- Stokes 
equations (|l|2|3p and have in particular obtained the coefficients of the CGL ([TP]) for realistic 
fluid parameters [22] . To this end we have performed a numerical linear stability analysis of the 
steady hexagons employing a Galerkin expansion of the velocity and temperature fields for the 
Floquet analysis. 

The linear coefficient a = a r + iai is readily obtained from the growth rate and frequency 
of the perturbation associated with the Hopf bifurcation. The calculation of the diffusion and 
dispersion coefficient d = d r + i di is a bit more involved. It captures the dependence of the 
complex growth rate a of the oscillatory Hopf mode on slow spatial modulations of its amplitude, 
a = —s 2 d, where s is the wavenumber of the modulation, i.e. the Floquet exponent. Note, that 
the dependence of the growth rate of the translation mode (rather than the Hopf mode) on the 
Floquet parameter would yield the phase diffusion coefficients. We choose a discrete number 
of modulation wavevectors s typically in the range 0.01 < \s\/q c < 1 for the computations and 
obtain d by a fit of a to the equation a = — s 2 d. FigJH shows the ratio di/d r for To = 14°C and 
h = 4.92 mm as a function of the deviation of wavenumber of the underlying hexagons from 
the critical wavenumber, Aq = q — q c . Within the coupled Ginzburg-Landau equations di/d r 
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Figure 4. (Color online). Dependence of the ratio di/d r on the wavenumber Aq = q — q c . To = 14° C, 
Q = 65, non-Boussinesq coefficients 7§ = 0.0009, 71 = 0.2937, 7§ = -0.1681, 73 = 0.0215, 74 = 
-0.0022, and Prandtl number Pr = 8.2. 



vanishes at the bandcenter, Aq = 0, and is quadratic in Aq. The fully nonlinear computations 
of the Navier-Stokes equations confirm this behavior quite well. 

For the calculation of the non-linear coefficient c = c r + zci we perform direct simulations 
of the Navier-Stokes equations (p]l4|) in a small box of size L = 2 x 2it/q c to avoid side-band 
perturbations and spatial modulations of the oscillation amplitude. Therefore, from the real 
part dH/dt = a r \TL\ + c r \H\ 3 we extract c r by fitting the time derivative of the amplitude of 
the Hopf mode as a function of its amplitude (notice that we extracted a r in the previous step 
through a stability linear analysis). Finally, from the imaginary part uo = ai + Ci\H\ 2 we obtain 
C{ by fitting the oscillatory frequency of the Hopf mode as a function of its amplitude. 

For To = 14° C and h = 4.92mm with £2 = 65 we obtain coefficients that match closely the 
expectations based on the weakly nonlinear theory provided by the three coupled Ginzburg- 
Landau equations. In particular, we confirm that the system should exhibit bist ability between 
homogeneous oscillations and defect chaos. Direct simulations in sufficiently large system con- 
firmed this prediction. When starting with ordered hexagons and small perturbations that break 
the chiral symmetry, i.e. the three dominant wavevectors making up the hexagon pattern are 
given slightly different amplitudes, we obtain stable homogeneous oscillations, i.e. all hexagons 
whirl in phase with each other. To obtain the defect-chaos, defects need to be introduced into 
the complex oscillation amplitude. Note, however, that preferably there should be no defects in 
the underlying hexagon pattern. We therefore start with initial conditions in which the complex 
amplitude exhibits phase winding and its phase changes by 2tt across the system. In our simu- 
lations we find that this phase-modulated state is linearly unstable to side-band perturbations 
of the complex oscillation amplitudes which lead to defects in that amplitude. Through their 
wavenumber selection the defects induce the creation of further defects, without any significant 
deformations of the hexagon lattice. 

For somewhat stronger non-Boussinesq effects (To = 12°C with 4.92mm) we find that 
the Hopf bifurcation becomes subcritical. Since within the three coupled Ginzburg-Landau 
equations the Hopf bifurcation is always supercritical, this shows that the Hopf bifurcation 
occurs now beyond their regime of validity. Nevertheless, near the parameter values for which 
the bifurcation changes direction the whirling hexagons can still be described by a complex 
Ginzburg-Landau equation if a quintic saturating term — g\H 4 l~l is retained. Following the 
same procedure as before we determine the coefficients of this quintic equation and obtain 
d = 22.4 + 0.389z, c = -12.96 + 84.7z, g = 42.53 + 17.20i. The dependence of di/d r on the 
hexagon wavenumber is shown in Figj5l In contrast to the case To = 14° C shown in Fig 21 the 
dependence is now essentially linear near q c , reflecting again the fact that the Hopf bifurcation 
occurs at larger amplitudes for which the coupled Ginzburg-Landau equations are no longer 
valid. The inset of Figj5] shows that the dependence of co with respect to TL 2 is almost linear, 
and thus the quintic contribution gi\TL\ A to the frequency is very small in this regime. 

Direct simulations of the Navier-Stokes equations reveal that in this regime typical states 
exhibit localized unsteady bursts in the oscillation amplitude as shown in FigEl The white lines 
give the hexagonal cells, showing the deformations due to the whirling, and the magnitude \TL\ 
of the oscillation amplitude is shown in the color scale (red: maximal, blue: zero). We extract 
the complex oscillation amplitude TL from the simulation by demodulation with respect to the 
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Figure 5. Dependence of the ratio di/d r on the wavenumber Aq = q — q c (main figure), and dependence 
of the frequency of the Hopf mode on its amplitude (inset). The ratio di/d r is extracted from linear 
stability calculations and the frequency from numerical simulations in a small box of size L = 2x2 tv /q c , 
with T = 12° C, Q = 65, non-Boussinesq coefficients 7§ = 0.00100, 7? = 0.4885, 7§ = -0.2281, 
7§ = 0.0287, 7! = -0.00320, and Prandtl number Pr = 8.7. 




Figure 6. (Color online). Localized bursting of the oscillation amplitude. Color gives magnitude of 
oscillation amplitude (red: maximal, blue: zero). White lines indicate underlying hexagon pattern and 
show the deformation of the hexagon cells by the whirling. 



wavenumber of the hexagons and the dominant frequency of the oscillation. To do so, we first 
evaluate the Fourier transform of the time series given by the amplitude at each point of the 
numerical grid, and then shift the resulting temporal Fourier spectrum by the Hopf frequency. 
In this way we capture the slow time evolution. To obtain the slow spatial modulations, we 
compute the spatial Fourier transform of the shifted temporal Fourier transform, shift it such 
that one of the six dominant peaks making up the hexagonal lattice is moved to the origin, and 
remove the spectral content beyond a radius of V3q c /2. 

To assess the validity of the quintic complex Ginzburg-Landau equation we simulate it 
numerically using the parameters obtained from the simulations of the Navier-Stokes equation. 
We find that using the same system size the quintic CGL does not exhibit the bursting found 
in the Navier-Stokes simulations since the homogeneous oscillations become unstable to side- 
band perturbations only for smaller wavenumbers. If, however, the system size is chosen larger, 
similar bursting is obtained, as shown in FigJTl A more careful analysis shows that in this regime 
the coupling to the phase equations for the deformations of the underlying hexagons cannot be 
ignored [22]. These deformations represent changes in the local wavenumber of the hexagons. 
As our numerical stability analysis of the steady hexagons shows, in this regime the growth rate 
of the Hopf mode increases significantly with decreasing hexagon wavenumber. Thus, due to the 
deformation of the hexagons by the whirling the local growth rate within a burst is enhanced 
and allows this burst to persist. 

Increasing the Rayleigh number, the steady hexagons again become stable, as in the non- 
rotating case discussed above. The corresponding bifurcation diagram is shown in FigJHl Nev- 
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Figure 7. (Color online). Localized bursting of the oscillation amplitude in the quintic complex 
Ginzburg-Landau equation. Parameters correspond to those of Fig|6] except for a larger system size. 
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Figure 8. (Color online). Bifurcation diagram of non-Boussinesq convection with rotation for the 
subcritical case with To = 12° C and £2 — 65. Blue line corresponds to the minimal oscillation amplitude 
of the period- doubled orbit, the maximal oscillation amplitude is marked by the red line. 




Figure 9. Bursting state in water for To 
simulations of Eqs. Jl](4]) in a cell of size L - 
corresponds to time t — 663 t v and right t = 



= 12° C, Q — 65 and e = 0.7. The snapshots come from 
16 x 27i/q c using random initial conditions. Left snapshot 
699 t v . 



ertheless, for suitable initial conditions we still find irregular bursting of the whirling activity 
for e = 0.7, but not for e = 0.9. Two snapshots depicting whirling localized within domains 
of varying sizes are shown in Fig|9l Compared to the bursting obtained for e = 0.5 [22], the 
fluctuations in the whirling activity, I(t) = N~ x Jj^| >0 5 |^| \TL\dxdy with N denoting the 
temporal mean of the integral, are much smaller (FigfTD)). Nevertheless, persistently shifting 
domains of whirling hexagons are still apparent in Fig|9l 



Increasing e further, we find a period-doubling bifurcation at e c± 1.3 as shown in FigJHl which 
gives the larger (red) and smaller (blue) peak amplitude of the oscillations. Fig {TH presents a 
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Figure 10. Temporal evolution of the whirling activity lit) — N 1 f in ,,^ n K]n ,\ 1^1 dxdy of the 
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bursting state (see Fig(9}. 
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Figure 11. Period doubling after a Hopf bifurcation in the amplitude of one of the three modes that 
forms a hexagonal lattice. The simulations are made in a small cell of size of L = 2 x 2n/q c with 
To = 12, e= 1.5, and Q = 65. 



time series of the amplitude of one of the Fourier modes that forms the hexagonal lattice for 
e = 1.5. The time series is obtained from the simulation of equations (1 — 4) in a small cell of size 
L = 2 x 2ir/q c . Whether there are any further period doubling events for larger e remains to be 
investigated. So far we have obtained this state only in simulations of small systems in which 
side-band perturbations of the oscillation mode as well as the hexagon pattern are negligible. 

As the non-Boussinesq effects are increased yet further [23] , the range over which the steady 
hexagons are unstable to the whirling becomes smaller and eventually the steady hexagons are 
linearly stable over the whole range of Rayleigh numbers considered, < e < 1.2. However, since 
the Hopf bifurcation is subcritical in this regime solutions with whirling hexagons still exist. In 
fact, the whirling activity, which does deform the underlying hexagon pattern |3Q|31|33j . is so 
strong that it breaks up the hexagonal lattice, introducing many defects. In weakly disordered 
hexagon patterns the dominant defects are typically penta-hepta defects which consist of two 
adjacent convection cells that have 5 and 7 immediate neighbors, respectively. The state result- 
ing from the strong whirling activity is, however, considerably more disordered and exhibits 
not only such cells with coordination numbers 5 and 7, but also a small number of cells with 
coordination numbers 4 and 8, respectively [23] . 

For intermediate Rayleigh numbers the spatio-temporally chaotic whirling activity arises 
only transiently. During that transient one can qualitatively discern a mutual reinforcement of 
the defect formation and the whirling activity. This is illustrated in the two snapshots shown in 
FigfT2l Whirling activity emanates from the defects in the hexagonal lattice (marked by a red 
and a blue circle). It in turn tends to induce other defects, see right panel. For the parameters 
chosen in Fig[l2]the whirling emanating from the defects is not strong enough and the various 
defects in the lattice eventually annihilate each other and the pattern becomes steady. If the 
mutual reinforcement is strong enough, however, persistent activity can arise. For e = 1 the 
whirling invades the steady state and in systems of size L = 16 x 2ir/q with q = 4.5 we find that 
the chaotic activity persists as long as we could simulate the system (more than 300 vertical 
thermal diffusion times, which corresponds to more than 5 hours physical time). 



Will be inserted by the editor 
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Figure 12. (Color online). Interplay between whirling activity and defect formation. Two snapshots 
for e = 0.87 and h = 4.6mm, To = 12° C £2 = 65. Penta-hepta defects (one is marked by blue and red 
circle) induce whirling and whirling tends to induce defects. 



5 Conclusion 



In this paper we have elaborated on some recent results for hexagon convection in non-Boussinesq 
convection [21(22(20] . In the non-rotating case we focused on the reentrance of hexagons as the 
Rayleigh number is increased. We contrasted the reentrance found in our computations for 
water with that found in experiments using a layer of SF6. The reentrant hexagons obtained 
through the mechanism elucidated with our computations generically are expected to inherit 
the low wavenumber that is characteristic of Boussinesq hexagons [16] . In contrast, the reen- 
trant hexagons found in the SF6 experiments have a wavenumber close to q c [15], suggesting 
that their reentrance is due to a different mechanism. 

For stronger non-Boussinesq effects the hexagons do not exhibit any linear instability to 
rolls (or the mixed mode) and, in fact, even defects arising through side-band instabilities 
of the hexagons do not necessarily lead to a transition to rolls. Instead the pattern can heal 
to a regular hexagon pattern with lower wavenumber. In convection in gases, which have a 
low Prandtl number, the side-band instabilities limiting the stable wavenumber band can be 
different at the \ow-q and the high-g side, which leads to transients exhibiting different degrees 
of disorder [20] . 

Similar looking reentrant hexagons have also been reported in a chemical model system, the 
two-dimensional Brusselator system, for which a series of numerical simulations were performed 
[29] . The initial Turing instability off the homogeneous state yields hexagon patterns. With in- 
creasing control parameter these hexagons lose their stability to stripe-like patterns. For yet 
larger values of the control parameter the hexagons regain their stability. In contrast to the 
behavior found in convection in this chemical model the hexagons near onset are up-hexagons 
(7r-hexagons) while the reentrant hexagons are down- hexagons (0-hexagons). These results are 
modeled with an amplitude equation whose quadratic coefficient depends on the control pa- 
rameter. Such a model has also has been used for hexagons in magnetic fluids [34] and in our 
previous calculations for water [2T] . 

For rotating systems we discussed the whirling hexagons that arise from a Hopf bifurcation 
that takes the place of the steady bifurcation to a mixed mode in the non-rotating case. In 
the weakly non-Boussinesq case this transition occurs for very small convection amplitudes 
and the Hopf bifurcation is supercritical, as predicted by computations based on three coupled 
Ginzburg-Landau equations describing weakly nonlinear hexagon patterns with rotation [30] . 
It turns out that to leading order the oscillations can be described by the two-dimensional 
complex Ginzburg-Landau equation. Moreover, the resulting oscillations (whirling hexagons) 
generically arise in the regime in which defect chaos persists despite the existence of stable 
plane wave and homogeneously oscillating solutions. For stronger non-Boussinesq effects the 
bifurcation becomes subcritical and we find localized bursting in the whirling amplitude and 
also whirling activity that is strong enough to tear the hexagonal lattice apart. 
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